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Abstract 

We provide accurate assessments of the consequences of violations of self-consistency in Hartree- 
Fock (HF) based random phase approximation (RPA) calculations of the centroid energy -Ecen of 
isoscalar and isovector giant resonances of multi-polarities L = — 3 in a wide range of nuclei. 
This is done by carrying out highly accurate HF-RPA calculations neglecting the particle-hole (ph) 
spin-orbit or Coulomb interaction in the RPA and comparing with the fully self-consistent HF-RPA 
results. We find that the shifts in the value of Ecen due to self-consistency violation associated with 
the spin-orbit and Coulomb interactions are comparable or larger than the current experimental 
errors in Ecen- 

PACS numbers: 24.30. Cz, 21.65. +f, 21.60.Jz 
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I. INTRODUCTION 



The study of collective modes in nuclei provide very important information for under- 
standing the structural and bulk properties of nuclear systems. For example, the isovector 
giant dipole (IVGDR) mode is sensitive to the symmetry energy, the centroid energy Ecen of 
the isoscalar giant monopole resonance (ISGMR) is related to the value of the incompress- 
ibility modulus K^o of symmetric nuclear matter , and low lying collective states give 
a clue on the nuclear shell structure. These quantities are important ingredients not only 
for the description of finite nuclei but also for the study of heavy-ion collisions, supernovae 
and neutron stars. Recent developments in high precision experimental facilities make it 
possible to measure the centroid energy of the ISGMR with an error of SEcen ~ 0.1 — 0.3 
MeV 0,0]. Using the approximate relation {6Koo)/Koo = '^{SEcen)/ Ecen and, for example, 
the recent experimental value of Ecen = 13.96 ± 0.20 MeV for the ISGMR in ^°^Pb, one has 
an error of 6Koo = 6 — 9 MeV for Koo = 200 — 300 MeV. This enhanced precision calls for a 
critical accuracy check at the side of the the theoretical calculations with the goal that the 
error in the calculated value of E^en is less than the experimental error. 

The basic theory for the microscopic description of different modes of giant resonances is 
the Hartree-Fock(HF) based random phase approximation (RPA) ^ 1^ . Although this ap- 
proach is conceptually well understood, actual calculations make compromises for reasons of 
simplicity or numerical expense. One can hardly avoid limitations in the numerical basis and 
RPA phase space. Furthermore, most of the presently available HF-RPA calculations omit 
the painful to evaluate pieces of the RPA residual interaction, namely its spin-orbit and/or 
Coulomb parts. We will call that the self-consistency violation (SCV) in the following. It 
is obvious that a very accurate calculation within HF-RPA demands a sufficiently complete 
basis and in particular self-consistency, i.e., using exactly the same pieces in the residual 
interaction that have been used in the underlying HF calculation. Unfortunately, apart from 
»™e M,y .e.f-c„„...„t calculate. flB&&0, e..t.g HF-RPA calculations 
are contaminated by SCV. First assessments on the effects of SCV on the strength functions 
of giant resonances were obtained in Refs. Q, Q, Q In Ref. 0, 

results of elaborate 

studies of the effects of SCV on the constrained energy (Econ) and scaling energy (Eg) have 
been reported only for the ISGMR. It is to be noted that the full self-consistent values of 
Econ and Eg were obtained using the constrained Hartree-Fock (CHF) approach. It was 



2 



pointed out that the SCV concerning spin-orbit and Coulomb interactions may cause 
an error in Econ of the ISGMR which becomes larger than 1 MeV, i.e. as large as 5 times 
the experimental error. This calls for a systematic assessment of the effects of SCV on the 
excitation energies of other modes as well. Recently in Ref. ^, fully self-consistent calcu- 
lations of strength function were reported and the shift of the peak of the ISGMR strength 
function due to SCV has been discussed for ^^°Sn. In Ref. the effect of SCV has been 
assessed for the IVGDR in ^^O using the fully self-consistent approach of small amplitude 
time dependent Hartree-Fock method. We will continue here with discussing a greater va- 
riety of modes, namely the sequence of isoscalar and isovector multi-pole resonances in the 
range L = — 3 and a few examples from low lying collective states. 

Before attacking the main task, it is worthwhile to put the various sources of uncertainties 
into perspective. The HF-RPA method optimizes the modes in the space of one-particle- 
one-hole (Iph) excitations. Correlations, i.e. 2ph and higher structures are not accounted for 
explicitlv. Such correlations have been very much discussed in the past, for reviews see e.g. 



)ncitiv. 

,0,0 



milallSl- The main effect is a coUisional broadening of the strength distributions. This can 



be accompanied by a certain shift of the resonance peak position which grows with excitation 
energy and can reach orders of 1 MeV for the rather high lying isovector modes (in the range 
above 20 MeV). However, the Skyrme forces employed in nuclear HP and RPA calculations 
are effective forces which incorporate already a great deal of correlations [2^. This reduces 
the correlation effects on the peak positions [lol I21I . We adopt the view that the net effect 
remains of order of a tenth of MeV for modes with moderate excitation energy around and 
below 15 MeV. A second crucial aspect concerns limitations within the HP-RPA approach 
itself. There are two major questions: the size of the RPA phase space and the handling of 
the particle continuum. We take care to use a sufficiently large expansion basis such that 
peak positions have converged to uncertainties below 0.1 MeV. The artificial discretization 
of the continuum has very little effect on the average peak positions but limits the spectral 
resolution with which one can reliably compute the strength distributions j^]. We use a 
large simulation box and properly adapted smoothing width. After all, there remain the 
effects of SCV to be studied. We demonstrate the accuracy of our fully self-consistent HP 
based RPA calculations of the strength functions of giant resonances by comparing (i) the 
RPA results with the corresponding ones of CHF for the case of the ISGMR and (ii) the 
total energy weighted strengths with the corresponding energy weighted sum rules (EWSR). 
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II. FORMALITIES 



We calculate the strength function, 

s{e) = J2e,\<o\f,\j>\' (1) 

j 

: oUowing the fully self-consistent method based on Q — P representation described in Ref. 
3, 3|- In Eq- ©; |0 > is the RPA ground state and the sum is over all RPA excited states 
\j > with the corresponding excitation energy Ej. The scattering operator is given by 

Fl = ^f{ri)YLoii), for isoscalar, (2) 

i 

Fl = —^f{rn))YLo{n) - — ^f{r{p))YLo{p), for isovector, (3) 

n p 

with /(r) = r^,r^ and for monopole, quadrupole and octopole, respectively. For the 
isovector dipole we have taken, /(r) = r, whereas for the isoscalar dipole we adopt the 
scattering operator /(r) = — | < > r to eliminate the contribution of spurious state 



mixing 



12| . Once we have the strength function, the energy moments can be calculated by. 



oo 



ruk = uj^S{u)duj. (4) 
Jo 

Then the centroid, constrained and scaling energies of the giant resonance are computed as. 



TP TP / A TP (K\ 

^cen = , J^con = \ , and Hjg = \ ■ [o) 



For the isoscalar (T = 0) operator in Eq. (j2)), the EWSR is given by 



m,(L,T^a) = ^^A{g,(r)), (6) 

with 

{giir)) = -^j gL{r)p{r)'iTTr^dr, (7) 
where p(r) is the HF ground state matter density distribution and 

3.M^(|)%i(i + i)(9'. (8) 

For the isovector (T = 1) operator of Eq. 0, we have, 

NZ 

mi(L, T = 1) = -^mi (L, T = 0) [1 + k - K^p] (9) 



where k is the enhancement factor due to the momentum dependence of the effective nucleon- 
nucleon interaction, and is given by, 

_ (1/2) + Xi/2) + t^jl + X2/2)] 2 / gLir)ppir)pr,ir)A7Tr'dr 

{hy2m){4:NZ/A^) J gLir)p{r)ATTr^dr ' ^ ^ 

where ti and Xi are the parameters of the Skyrme interaction. The correction K^p, which 
arises because of the difference in the profiles of the neutrons and protons density distribution 
( ie., since p„(r) - pp{r) ^ ^^p(r)), is given by, 

_ {N-Z) A JgL{r)[Zpn{r)-Npp{r)]47rr^dr 
""^ A NZ J gL{r)p{r)47rr^dr ' ^ ^ 

The fully self-consistent moment m_i for the ISGMR can also be calculated using the 
constrained Hartree-Fock method Q], 



1 Id < > 



rri-i = — -- 



1 Id^E^ 



An 2 dX^ 



(12) 

A=0 



An 2 d\ 

where the mean square radius < >a and the energies E\ are obtained using the solution 
of the constrained Hamiltonian H\ = H — Xr^. The self-consistent energy moment mi can 
be evaluated |3 using the EWSR of Eqs. (0)-® with / = r^. 



1 h"^ 

mi = AA < >, (13) 

47r2m ' ^ ^ 

where m is the mass of the nucleon, A is the mass number and < > is the mean square 
radius calculated using the ground state HF wave function. Thus, using Eqs. flT^ and (O, 
one can get the self-consistent value of Econ which can also be used to check the accuracy of 
the RPA calculations (in particular, the completeness of the RPA phase space) . For plotting 
purpose, we employ a Lorentzian smearing of the strength function of Eq. ^ obtaining, 

where F is the smearing parameter taken to be 2 MeV or larger. Note, however, that for the 
evaluations of the energy moments, Eq. and the energies of various giant resonances, 
we use a very small value for F (< 0.2 MeV). 

Since we shall be investigating the magnitude of errors in the calculation of the observ- 
ables for giant resonances due to the violation of the self-consistency, it is necessary to pay 
attention to the accuracy of HF and RPA calculations. We have taken a box of size 30 fm. 
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This large computational box allows a spectral resolution of the strength functions in the 



particle continuum of about 0.8 MeV which is by far sufficient for our present inves- 
tigation. The mesh size is taken to be 0.3 fm. The flexible formulation of the RPA [3] in 
the Q-P representation allows one to include Iph excitations up to very high energy, just 
jy taking the detailed Iph space only for transitions up to few major single particle shell 
15[. We have taken the detailed Iph space up to 4 particle major shell and checked the 
convergence of the results by comparing with those obtained with 6 particle major shell. We 
have adopted the SGII jl^ interaction for the HF based RPA calculations and the Slater 
approximation has been used for the Coulomb exchange term, consistently both in the HF 
and RPA calculations. For the CHF calculations, we have taken a box of 15 fm with mesh 
size 0.1 fm and calculated each derivative in Eq. (jl2p using flve point formula with the 
increment 0.02 in the constraining parameter A. 



III. RESULTS AND DISCUSSION 



In order to check the accuracy of our CHF calculations for m_i in ISGMR, we have given 
the values for ^°Ca and ^osp^ Table HI calculated from Eq. ()12p using the mean square 
radii, m„i(r), and the energy, m„i(e), methods, for different values of the constraining 
parameter A. For comparison, the values of m_i obtained from RPA calculations Eqs. 
and (jH) are given in the last column. It can be seen that m_i(r) differs from m_i(e) by 0.7% 
in the case of A = 0.01 for ^°Ca. For A = 0.02, the deviations are very small (0.04%) for 
both the nuclei. In the following, we will use the value of A = 0.02 in our calculations. We 
add that for A = 0.02 the RPA values of m_i, obtained from Eq. using the integration 
range — 60 MeV, compare very well with the CHF values (within 0.6%). This clearly 
demonstrates the high accuracy of our CHF and HF-RPA calculations. 

A very important necessary condition for a fully self-consistent HF based RPA calcula- 
tions for giant resonances is to obtain the full EWSR from the calculated RPA strength 
function. In order to see how this necessary condition is fulflUed in our present calculation, 
we compare in the Table ITTl the mi values calculated using the RPA strength functions in 
Eq. (jH, {mi{RPA)) with the corresponding EWSR obtained from Eqs. ^ and (jH)) for the 
nuclei ^°Ca, ^°Zr and ^°^Pb for L = — 3 and T = and 1. It is seen that for all of these 
nuclei and for all the modes, the deviations of mi {RPA) from the corresponding EWSR are 



6 



very small (less than 0.3%). This once again shows the high accuracy of our HF-RPA calcu- 
lations. We note that for the SGII interaction we have n = 0.314, 0.381, 0.314 and 0.253 and 
the values of n^p = 0.010, 0.000, 0.010 and 0.024 for the isovector L = 0, 1, 2 and 3 in ^ospb, 
respectively. We point out that the correction term in O for isovector modes which is 
usually missing in the literature is not negligible for asymmetric nuclei. As an example, for 
^°^Pb, L = 3, Knp has an effect of 2% in the calculation of EWSR. The effect of Unp will 
be more significant for nuclei near the drip lines because of the large difference between the 
neutron and proton density distributions and also for the large asymmetry (A^ — Z)/A (see 
Eq. (HH)). 

In Table IllH we have given the energies Econ, Eg and Ecen of the ISGMR for a host 
of nuclei ranging from very light to heavy, including some neutron or proton rich nuclei. 
The second column of Table IIIII indicates the kind of calculation — SC, LS or CO. The 
symbol SC corresponds to the fully self-consistent calculation, i.e., both the spin-orbit and 
Coulomb interactions are present in HF as well as in RPA. The spin-orbit violation (LS) 
means that the spin-orbit and Coulomb interactions are present in HF but the ph spin-orbit 
interaction is missing in the RPA calculation. Finally, the Coulomb violation CO means that 
the spin-orbit and Coulomb interactions are present in HF but the ph Coulomb interaction 
is dropped out from the RPA calculations. The constrained energies Econ calculated with 
the CHF approach using Eqs. (fT^ and (fTSj) are presented in the third column. The values 
of Econi Es and Ecen of the ISGMR calculated in HF based RPA, are given in the next three 
columns. The range of integration for calculating the ISGMR energies from the strength 
function (Eqs. (jH) and ©) is 0-60 MeV. Comparison of Econ obtained from the fully self- 
consistent (SC) HF-RPA calculation of the strength function (forth column) with those 
obtained from the constrained HF method (third column) shows a maximum difference of 
0.06 MeV. This indicates that the accuracy of our RPA calculation is very good. It is to 
be noted that the effect of the SCV of LS and CO are similar in Econ and Eccn-, but a little 
different in Es. The shift of the centroid energy 5 Ecen due to LS or CO violation are listed 
in the last column. The effects of violation for the spin-orbit (LS) interaction are very small 
for ^^O, ^°'^°Ca and 80,ii02j~ i^^^ robust in other nuclei, particularly for ^""Sn (~1.2 MeV) 
and ^^Ni (~1.8 MeV). This can be understood because the effect of the omission of the 
ph spin-orbit interaction is expected to show up prominently only for the spin unsaturated 
nuclei like ^^Ni and ^"''Sn here. Hence, for the spin closed ^^O, ^°'^°Ca and 80,ii02j^ nuclei. 
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the effects of violation of spin-orbit interaction are very small. On the other hand, the effect 
of the CO violation seems to be dependent on the position of the nuclei in the Segre chart. 

The energies of the ISGMR are obtained from the ratios of the moments in Eq. (jSj). In 
Table HVl we give for ^°Ca and ^osp^ the values of mi, mo and m_i obtained from the fully 
self-consistent HF-RPA calculations (SC) along with those values obtained from calculations 
where self-consistency is violated due to the omission of ph spin-orbit (LS) and ph Coulomb 
(CO) interactions. Comparison of the mi values calculated in fully self-consistent RPA with 
the EWSR obtained from Eq. p3p (given in column 3) shows that the RPA values are 
smaller than the corresponding EWSR values by less than 0.3%. The violations of self- 
consistency (LS or CO) has little effect on values of mi (0.04%) while the effects of SCV 
are considerably larger in m_i and mo which mainly cause the shift in the giant resonance 
energies. This shows that obtaining a value of mi in a RPA calculation which is very close 
to the EWSR does not indicate an accurate and fully self-consistent implementation of HF 
based RPA theory (see also jl^). 

In theoretical study of giant resonances, the evaluation of the strength function S{E) 
is needed to compute the energy moments (jH) and hence the centroid energy, constrained 
energy, etc. Therefore, it is worth seeing how the strength function as such gets affected by 
the self consistency violation. In Fig. 1, we display the variation of the S{E) with energy 
E for isoscalar excitations of different multi-polarities (L = — 3), for the nucleus ^osp]^ 
a representative case. We have smoothened all the strength functions by a Lorentzian with 
a smearing parameter F = 2 MeV (see Eq. (fT^ ) in order to wipe out artefacts from the 
discretization of the continuum and to provide a smooth curve for better comparison of the 
effects of SCV in S{E). In the top panel of Fig. 1 we have three curves for the ISGMR 
(L=0,T=0) — the fully self-consistent (SC) (solid line) result, with violation for ph spin-orbit 
interaction in RPA (LS) (dashed line) and the result with violation of ph Coulomb interaction 
in RPA (CO) (open circle). All these three curves have similar single-peaked structure. But 
it is significant to notice that the violation of LS has pushed the peak to a higher energy 
by almost by 0.7 MeV, whereas the position of the peak of the curve has moved toward the 
lower energy side by around 0.4 MeV due to the self-consistency violation for the Coulomb 
interaction. This implies that the self consistency violation may cause an uncertainty in the 
calculated value of the centroid energy Ecen of ISGMR for ^osp]^ which is twice in magnitude 
compared to the experimental one . This uncertainty, as we have discussed above, would 
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impose an error of about 20 MeV in the prediction of the value of K^o- Since the effects 
of LS and CO violations are opposite to each other, there would be a partial cancellation 
if one does not take both of the ph LS and CO interactions simultaneously in the RPA 
calculations. For the isoscalar giant dipole (L=1,T=0) resonance (ISGDR), the effects of 
neglecting the Coulomb and spin-orbit terms in the ph interaction are similar to those for 
ISGMR. Violation of LS leads to a push of the peak of S{E) toward a higher energy whereas 
the violation of CO pulls the peak to a lower energy. But the magnitude of the shift of 
the peak is much smaller in comparison with the ISGMR. Specifically, the ph Coulomb 
interaction shows very little effect. We display similar curves for the higher multi-polarities, 
L=2 and 3, in the lower two panels. A shift of the position of the peak of S{E) is observed 
once we neglect either the ph Coulomb or spin-orbit interactions in the RPA calculation. 
In both cases of the quadrupole and octopole giant resonances, the effect of the violation of 
self consistency acts in the same direction which is in contrary to the ISGMR and ISGDR 
cases. Therefore, the total shift in the peak energy of S{E) will be very significant if one 
omits both Coulomb and spin-orbit ph interactions in the RPA calculations. Note also the 
significant shifts in the energies and strengths of the low lying isoscalar L = 1 — 3 states due 
to SCV. 

In Fig. 2, we plot the strength functions for the isovector modes for ^osp]^ the same way 
as we have done in Fig. 1 for the isoscalar mode. Here we have used in the calculations the 
same parameters as those taken for the isoscalar modes except for the smearing parameter 
r = 10 MeV for the isovector giant monopole resonance (IVGMR). It is clear from the figure 
that in the isovector modes (for all considered multi-polarities) the strength functions are 
almost insensitive to the self consistency violation due to the omission of the ph spin-orbit 
interaction in the RPA calculation. On the other hand, the absence of the ph Coulomb 
interaction in the RPA calculation pushes the strength function significantly toward lower 
energy and this is most prominent in the isovector dipole mode where the shift is around 
0.4 MeV. 

In Table IVj we present, for a wide range of nuclei, the self consistent centroid energies 
Ecen{SC) and their shifts 5Els and SEco due to the self-consistency violations for ph spin- 
orbit and Coulomb interactions in the RPA calculations, respectively. We have given results 
for different multi-polarities (L = 0,1,2 and 3) for both the isoscalar and the isovector modes. 
It is evident from the table that for isoscalar modes, the effect of LS or CO violations are most 
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prominent for monopole resonance which are almost 2 to 3 times larger than those for the 
other multi-polarities. Note that if we drop the ph LS and CO interactions simultaneously in 
the RPA calculations for ^°^Pb ISGMR, the shift of the centroid energy is SEls.co = —0.30 
MeV which is comparable to the experimental uncertainty. The effects of the LS and CO 
self-consistency violations are found to be somewhat smaller for L = 2 and 3 in comparison 
with the ISGMR for all of these nuclei. But for the L = 2 and 3 modes, SEls and 6Eco 
are found to be of the same sign. Therefore, their combined effect on the centroid energy 
is significant. The SCV associated with the Coloumb interaction (in the RPA level) affects 
Ecen considerably in the dipole and monopole modes of the isovector channel. 

Finally, we have checked the effect of SCV on low lying (isoscalar) collective states. Test 
cases were the lowest 3~ in ^^Zi and the lowest 2+ in ^°^Pb. The low lying states are much 
more sensitive to the size of the RPA phase space. We have used here an expansion basis 
up to about 2000 MeV to ensure sufficient convergence. The results are shown in Table IVIl 
The effects are of the same order as found for the giant resonances. However, for these low 
lying states a mismatch of about 0.5 MeV is large compared to the total excitation energy. 
A fully consistent calculation is compulsory for studying spectra in that energy range. 

IV. CONCLUSION 

In summary, we have carried out highly accurate fully self-consistent Hartree-Fock (HF) 
based random phase approximation (RPA) calculations for the strength functions of isoscalar 
and isovector L = — 3 modes in a wide range of nuclei. We have quantified very accurately 
the effects of self-consistency violations in the calculations of the energies of giant resonances 
of nuclei within the HF based RPA. We have studied the cases of SCV due to the omission 
of the spin-orbit (LS) or/and Coulomb (CO) ph interactions and mainly focus on their 
effects on the centroid energy Ecen- Here we consider both isoscalar and isovector modes of 
multi-polarities L = — 3. It is found, for the wide range of nuclei considered here, that 
the effects of violations of self-consistency due to the ph LS or CO interactions are most 
significant for the ISGMR. For the ISGMR, the absence of the ph LS interaction tends to 
increase Ecen, whereas the violation due to ph CO interaction decreases Ecen- For the spin 
unsaturated nuclei (such as ^^Ni and ^''^Sn), the shift in Ecen is robust (~ L5 MeV) which is 
almost 5 times larger than the experimental uncertainty. For other higher multi-polarities. 
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the individual effect of the ph LS and CO interactions are somewhat smaller than those for 
the ISGMR. But for the quadrupole and octopole modes, the LS and CO self-consistency 
violations both tends to reduce the centroid energy. Hence, the effect of SCV on Ecen in 
these modes are significant (0.3 — 0.6 MeV) if one neglects the ph spin-orbit and Coulomb 
interactions simultaneously in the RPA calculation. 
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TABLE I: The values of m_i of Ca and Pb for ISGMR obtained from the mean square radii, 
m_i(r), and energy, ?n_i(e) methods (Eq. ()12() ). are presented for different values of the increment 
of the constraining parameter A. For comparison, the RPA values are given in the last column. 



Nucleus 




CHF 




RPA 




A 


m_i(r) 


m_i(e) 






(fm-2 MeV) 


(fm^ MeV-i) 


(fm" MeV-i) 


(fm'^ MeV-i) 


40Ca 


0.01 


6.512 


6.467 






0.02 


6.469 


6.466 


6.426 




0.03 


6.466 


6.463 




208pb 


0.01 


230.01 


230.57 






0.02 


230.57 


230.50 


230.83 




0.03 


230.27 


229.89 





TABLE II: Comparison of mi calculated from RPA strength function S(E), mi{RPA) (Eq. Q 
integrated up to Emax ~ 100 MeV), with those obtained from the energy weighted sum rules 
(EWSR) in Eqs. © and ®. The ratio R = mi{RPA)/EWSR indicates that the maximum 
deviation of mi{RPA) from EWSR in this table is less than 0.3%. 



Mode 


Ca40 


Zr90 


Pb208 


L 


T 


EWSR 


R 


EWSR 


R 


EWSR R 








2889 


0.9983 


10505 


0.9994 


41872 0.9971 




1 


896.3 


0.9992 


3330 


0.9994 


13041 0.9998 


1 





57253 


0.9990 


289907 


0.9995 


1823110 0.9999 




1 


64.62 


0.9999 


148.9 


1.0001 


337.8 0.9999 


2 





7222 


1.0001 


26262 


1.0007 


104681 1.0015 




1 


2241 


1.0000 


8326 


0.9998 


32604 1.0002 


3 





238240 0.9996 


1300645 0.9999 


8584813 1.0003 




1 


69328 


0.9994 


389266 


0.9996 


2519881 1.0001 
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TABLE III: The constrained {Econ), scaling (Es) and the centroid energies (Ecen), in MeV, of the 
ISGMR of various nuclei are given for fully self-consistent HF-RPA calculations (SC) along with 
those obtained for the calculations where the self-consistency is violated due to the neglecting of the 
ph spin-orbit (LS) and Coulomb (CO) interactions in RPA. We adopt the range {u)i—uJ2) = (0 — 60) 
MeV for the energy integration. The effects of the SCV on the centroid energies are given in the 
last column. We have also given in Column 3, the constrained energies obtained from fully self- 
consistent constrained Hartree-Fock (CHF) calculations. The Skyrme interaction SGII was 
used. 



Nucleus 


9y 


CHF 




RPA 








Econ 


Econ 


E, 


E 


^Ecen 


160 


SC 


24.025 


23.985 


26.846 


24.584 






LS 




23.710 


26.905 


24.384 


0.200 




CO 




23.890 


26.722 


24.481 


0.103 


40 Ca 


SC 


21.131 


21.186 


22.516 


21.469 






LS 




21.026 


22.605 


21.362 


0.107 




CO 




20.960 


22.258 


21.234 


0.235 


60Ca 


SC 


16.203 


16.154 


18.804 


16.948 






LS 




15.956 


18.886 


16.812 


0.136 




CO 




16.093 


18.660 


16.865 


0.083 


56 Ni 


SC 


20.177 


20.179 


20.864 


20.311 






LS 




21.853 


23.056 


22.075 


-1.764 




CO 




19.883 


20.550 


20.009 


0.302 


80 Zr 


SC 


17.841 


17.873 


18.629 


18.050 






LS 




17.842 


18.752 


18.051 


-0.001 




CO 




17.527 


18.267 


17.699 


0.351 


90 Zr 


SC 


17.814 


17.832 


18.261 


17.914 






LS 




18.512 


19.233 


18.645 


-0.731 




CO 




17.532 


17.944 


17.609 


0.305 


ii02r 


SC 


15.294 


15.248 


16.410 


15.585 






LS 




15.207 


16.596 


15.584 


0.001 




CO 




15.080 


16.160 


15.392 


0.193 


loogn 


SC 


17.156 


17.207 


17.645 


17.282 






LS 




18.351 


19.086 


18.488 


-1.206 




CO 




16.730 


17.189 


16.804 


0.478 


iiesn 


SC 


16.261 


16.288 


16.832 


16.400 






LS 




16.868 


17.659 


17.027 


-0.627 




CO 




15.961 


16.478 


16.062 


0.338 


i44Sm 


SC 


15.276 


15.305 


15.806 


15.372 






LS 




15.991 


16.740 


16.100 


-0.728 




CO 




14.919 


15.459 


14.987 


0.385 


208pb 


SC 


13.475 


13.454 


13.910 


13.522 






LS 




14.080 


14.767 


14.196 


-0.674 




CO 




13.079 


13.545 


13.142 


0.380 
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TABLE IV: The values of nik of Ca and Pb for ISGMR obtained from fully self-consistent 
RPA calculations (SC) along with those obtained from the calculations where the self-consistency 
is violated due to the neglecting of the ph spin-orbit (LS) and Coulomb (CO) interactions in RPA. 
The EWSR values of mi obtained from Eq. 1)131) are given in column 3. 



Nucleus 


9y 


EWSR 




RPA 








mi(Eq. (0) 


mi 




mo 






(fm^^ MeV-i) 


(fm^^ MeV) (fm^ MeV"^) 


(fm^) 


40Ca 


SC 


2888.7 


2884.4 


6.426 


134.35 




LS 




2884.6 


6.525 


135.03 




CO 




2885.9 


6.569 


135.91 


208pb 


SC 


41872.5 


41783.7 


230.83 


3090.03 




LS 




41782.6 


210.76 


2943.37 




CO 




41819.6 


244.49 


3182.06 
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TABLE V: Fully self-consistent HF-RPA results of centroid energies {Ecen{SC)) and their shifts (in 
MeV) due to self-consistency violations via spin-orbit {5Els = Ecen{SC) — Ecen{LS)) and Coulomb 
{6Eco = Ecen{SC) — Ecen{CO)) pli interactions are presented for different modes (L = — 3, 
T = 0, 1) of some nuclei of experimental interest. The ranges of energies of integrations {uji — UJ2) 
for the calculations of the centroid energies are given in columns 3 and 7. The Skyrme interaction 

sGii y was used. 



Nucleus 


L 




Isoscalar 








Isovector 












Ecen{SC) 


5Els 


5Eco 


UJl - UJ2 


Ecen{SC) 


6Els 


SEco 


40 Ca 





5-40 


21.31 


0.12 


0.21 


10-60 


32.87 


0.05 


0.17 




1 


22-60 


32.60 


-0.09 


0.08 


5-40 


19.03 


-0.01 


0.18 




2 


8-26 


17.16 


0.24 


0.07 


10-60 


30.46 


-0.01 


0.07 




3 


22-45 


31.48 


0.30 


0.02 


24-65 


39.87 


0.02 


0.02 


56 Ni 





5-40 


20.20 


-1.68 


0.30 


10-60 


33.74 


-0.29 


0.18 




1 


22-60 


32.78 


-0.54 


0.10 


5-40 


19.01 


-0.16 


0.21 




2 


8-26 


17.52 


0.32 


0.12 


10-60 


30.68 


-0.02 


0.10 




3 


22-45 


31.17 


0.55 


0.03 


24-65 


40.40 


0.14 


0.04 


90 Zr 





5-35 


17.89 


-0.68 


0.30 


10-55 


31.77 


-0.10 


0.27 




1 


18-50 


28.88 


-0.33 


0.09 


5-35 


16.78 


-0.03 


0.29 




2 


8-20 


14.64 


0.15 


0.17 


10-50 


27.53 


-0.05 


0.14 




3 


20-40 


27.31 


0.29 


0.05 


22-60 


36.52 


-0.01 


0.05 







5-35 


16.38 


-0.59 


0.34 


10-50 


30.26 


-0.20 


0.30 




1 


18-45 


27.39 


-0.33 


0.10 


5-35 


15.75 


-0.08 


0.32 




2 


8-18 


13.60 


0.14 


0.19 


10-45 


25.82 


-0.17 


0.16 




3 


18-32 


25.37 


0.40 


0.05 


22-55 


34.87 


-0.10 


0.05 


i44Sm 





5-35 


15.34 


-0.69 


0.39 


10-50 


29.92 


-0.08 


0.37 




1 


18-45 


26.42 


-0.33 


0.13 


5-35 


15.02 


-0.03 


0.39 




2 


8-18 


12.97 


0.11 


0.28 


10-45 


24.75 


-0.06 


0.22 




3 


18-32 


23.95 


0.14 


0.04 


20-50 


33.23 


-0.04 


0.10 


208pb 





5-30 


13.50 


-0.66 


0.38 


10-45 


27.79 


-0.08 


0.39 




1 


16-40 


24.04 


-0.35 


0.12 


5-30 


13.66 


-0.02 


0.43 




2 


8-16 


11.64 


0.11 


0.19 


10-40 


22.64 


-0.05 


0.25 




3 


15-28 


21.13 


0.24 


0.04 


20-50 


30.92 


-0.07 


0.08 



16 



TABLE VI: The effect of SCV for low lying isoscalar states in two selected cases. The Skyrme 
interaction SGII [1^1 was used. 



Nucleus 


L 


energies [MeV] 
SC LS CO LS-CO 


90 Zr 


3 


2.10 2.57 1.95 2.47 


208pb 


2 


4.48 5.06 4.38 4.99 
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FIG. 1: HF-RPA results for isoscalar strength functions of Pb for L = — 3 multi-polarities are 
displayed. SC (full line) corresponds to the fully self-consistent calculation where LS (dashed line) 
and CO (open circle) represent the calculations without the ph spin-orbit and Coulomb interactions 
in the RPA, respectively. A magnified giant resonance is shown in the inset (lowest panel) for L = 3. 
The Skyrme interaction SGII [iq] was used. 
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FIG. 2: Same as Fig^but for isovector strength functions. 
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